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ABSTRACT 

We present new mid-infrared images of the central region of the Orion Nebula using 
the newly commissioned SOFIA airborne telescope and its 5-40 /um camera FORCAST. 
The 37.1 /um images represent the highest resolution observations (< 4") ever obtained 
of this region at these wavelengths. After BN/KL (which is described in a separate letter 
in this issue), the dominant source at all wavelengths except 37.1 ^m is the Ney- Allen 
Nebula, a crescent-shaped extended source associated with 6^ D Ori. The morphology 
of the Ney- Allen nebula in our images is consistent with the interpretati on that it is 



ambient dust swept up by the stellar wind from 6^ D Ori, as suggested bv lSmith et al 



(|2005l ). Our observations also reveal emission from two "proplyds" (proto-planetary 
disks), and a few embedded young stellar objects (YSOs; IRc 9, and OMCIS IRSl, 
2, and 10). The spectral e nergy distribution for IRc 9 is presented and fitted with 



standard YSO models from iRobitaille et al.l (|2007l ) to constrain the total luminosity 



disk size, and envelope size. The diffuse, nebular emission we observe at all FORCAST 
wavelengths is most likely from the background photodissociation region (PDR) and 
shows structure that coincides roughly with Ha and [N II] emission. We conclude that 
the spatial variations in the diffuse emission are likely due to undulations in the surface 
of the background PDR. 

Subject headings: Infrared: ISM — Infrared: stars — ISM: individual objects (Orion 
Nebula[M42=NGC1976]) 



1. Introduction 

Because of its relative proximity and its favorable location on the near side of its parent 
molecular cloud, the Orion Nebula (M42, NGC1976) has long been studied as the archetype of 
massive star formation. Its infrared characteristics have been key to our understanding of the flow 
of energy and the role of dust in such a region. Extended thermal emission from warm dust in 
the vicinity of the Trapezium OB star cluster that powers the Orion Nebula was discovered at 



11.6 and 20 /xm by iNey &: AllenI (119691 ). More detailed observations of this "Ney- Allen Nebula" 
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were subsequently carried out between 8 and 13 /xm bv lGehrz et alj (|l975l ). who showed that the 
emission peaks on O"" D Ori (HD37023, B0.5V) and argued that the primary heating mechanism for 
the grains in this region is re sonantly trappe d Lyman-a radiation. Using somewhat higher spatial 
resolution observations (3"), iHavwardl (| 19941 ) observed the Orion Nebula at multiple wavelengths 
from 7.6 to 24.3 //m, and used the spectral energy d istributions (SEDs) on representative sightlines 
to confirm an earlier claim by iForrest et al.l (119761 ) of the presence of a second, cooler, extended 
emission component. Hayward showed that the cool component (110 K vs. 300 K for the warm 
component) peaks to the West of 9^ D Ori, and concluded that it arises from the back wall of the 
H II region blister and is primarily heated by radiation from 9^ C Ori (HD37022, 04-6pv). An 
overview of t he geometry of the Orion Nebula and the "Main Ionization Front" (MIF) is given by 
O'Delll (|200lh . 



Re cent high-resolution mid-IR observations of the Orion Nebul a have revealed a wealth of new 
detail (jSmith et al.l l2005l : iRobberto et al.l l2005l : iKassis et al.l 120061 ). The warm dust component 
is resolved largely into limb-brightened dust arcs around proplyds, and the brightest dust ar 
that a round 9^ T) Ori — i s either a consequence of that star being loca ted close to the backg round 
cloud ( Smith et al.ll2005l ) or is a photoevaporated circumstellar disk ( Robberto et al.ll2005l ). 



The 

remaining, mostly cool dust emission emerges as a complex of streaks, gaps, arcs, and protostellar 
jets. 

In this letter we report on a study of the Orion Nebula with the Stratospheric Observatory 
for Infrared Astronomy (SOFIA) using the mid-IR facility camera FORCAST. In this study, we 
extend the relatively high-resolution view of the diffuse emission in this region to longer wavelengths 
than have previously been reported, and focus on the Ney-Allen Nebula and the proplyds. Our 
report also includes brief descriptions of emission from the infrared source IRc 9 and the cluster 
of young stellar objects in OMCl South. Results for the active star formation region surrounding 
the B ecklin-Neugebauer sou rce and the Kleinmann-Low Nebula are reported separately in this 
issue (jPe Buizer et al.ll2012l ). 



2. Observations and Analysis 



The SOFIA observations of Orion were performed on the n ights of 30 November and 3 Decem- 
ber 2010 using the FORCAST instrument (| Adams et ahlboiol ). FORCAST is a dual-array mid- 
infrared camera capable of taking images with two filters simultaneously. The short wavelength 
camera utilizes a 256x256 pixel Si: As array optimized for 5-25 ^t m, while the long wavelength 
camera utilizes a 256x256 pixel Si:Sb array optimized for 25-40 //m (jHerter et al.ll2012l ). After cor- 
rection for focal plane distortion, FORCAST effectively samples at 0.768" per pixel, which yields 



^ Usage of the term "Ney-Allen Nebula" h as evolved to now re fer uniquely to this bright feature, which lies at the 
peak of the extended structure discovered bv lNev fc Allen l| 19691 ). 
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an approximately 3.4' x 3.2' instantaneous field-of-view. 

Observations were obtained in the 6.6, 7.7, 19.7, 31.5, and 37.1 filters at an aircraft altitude 
of 43000 ft. The chopping secondary of SOFIA was driven at 2 Hz, with a chop throw of 7'. Such 
large throws were necessary to assure that the sky reference beam was completely off the extended 
Orion nebula complex and sampling "empty" sky areas for background subtraction. For the same 
reasons the telescope was nodded 10' away from the Orion nebula. With this setup, Orion was 
sampled in only one of the four chop-nod beams, and therefore the final effective on-source exposure 
times for the observations were 150 s at 6.6 and 7.7 fim, 450 s at 19.7 iim, 300 s at 31.5 fim, and 
420 s at 37.1 fim. 

Because all FORCAST data were taken in the dichroic mode, we were able to determine 
relative astrometry using images from the two wavelengths that were observed simultaneously. Data 
were taken with the following filter pairs: 19.7 /im/37.1 fj,m, 19.7 /im/31.1 fim, 6.6 /im/37.1 fim, 
and 7.7 /xm/37.1 fim. Therefore the relative astrometry of all filters was determined precisely by 
bootstrapping from the 19.7 //m/37.1 //m pair. We estimate that the relative astrometry between 
filters is known to better than 0.5 pixels (~0.38"). Absolute astrometry for each image was fixed 
by assuming tha t the p eak of the BN source at 19.7 //m is the same as the peak at 7 mm given by 



Rodriguez et al.l ([20051), R.A.(J2000)=05h35ml4.119s, Dec(J2000)=-05d22m22.785s. 



The FORCAST images were calibrated to a flux density per pixel using a standard instrument 
response derived from measureme nts of standard sta rs and solar system objects over several short 



science and basic science flights (jHerter et al.ll2012l ). The color correction from a flat spectrum 



source to that of the Orion sources is < 2%. The estimated 3a uncertainty in the calibration due to 
variations in flat field, water vapor burden, and altitude is approximately ±20%. Flux calibrated 
images for the 7.7, 19.7, and 37.1 //m data are shown in Figs. [H [21 and[3l the images at 6.6 fxm 
suffered from poor image quality and are not shown. 

Our FORCAST images at 7.7 /im reveal a number of compact sources which are consistent 
with IRAC images of the Orion Nebula available from the Spitzer Legacy Archive and coincide with 
a number of known objects in the region (see Fig. [T|). Of the known sources detected at 7.7 /im, 
only LVl and IRc 9 are visible at 19.7 /im, and no known compact sources are detected at 31.5 and 
37.1 /im. Source fluxes (or upper limits) measured from our FORCAST images at 19.7, 31.5, and 
37.1 /im for IRc 9 and LVl are included in the SEDs shown in Figure HI 

A temperature map of the region was generated from the 19.7 and 37.1 //m images using a 
least-squares procedure to fit the observed surface brightnesses in each pixel to the model 

Iu = il- e-"-/("/"-/)) X B^iT,) (1) 

where B is the blackbody function. We have adopted the Mathis extinction law for the wavelength 
dependence of r (= r/r^e/) and r^e/ is evaluated for 19.7 /tm (Fig. [5]). The surface brightnesses at 
each wavelength were given equal weight. The model corresponds to the solution of the radiative 
transfer equation for emission/absorption of material along the line of sight. The median dust tem- 
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Orion Nebula: 7.7 ju.m 
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Fig. 1. — FORCAST images of the Orion Nebula at 7.7 fxm, including prominent known sources. 
Details for the BN/KL region (satur ated in this image) at the upper right are presented in a 
separate paper (jPe Buizer et al.ll2012l ). White crosses mark the locations of the Trapezium stars 
(9^ Ori). 




Fig. 2.— FORCAST image of the Orion Nebula at 19.7 fim. 
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Orion Nebula: 37.1 fj,m 




Fig. 3. — FORCAST image of the Orion Nebula at 37.1 /xm. The "Arc" of emission SW of the 
Trapezium and the bright feature NW of 9^ E Ori (dashed oval) are discussed in Section 13. 4[ 
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Fig. 4.— SEDs for IRc9 and LVl. Left Panel— SEP for IRc9 including our mid-IR F ORCAST 
meas urements, near-IR photometry from iMuench et alj ( 20021) . Spitzer IRAC fluxes ( Megeath 



2OI2I ) , and additional mid-IR ground-based measurements ([Smith et al.1 12005| ) . Photometry for 



IRc9 at 19.7 and 31.5 fim was d one using a 5. 4" apert ure; note that IRc9 was not observed at 
37.1 /xm. The best model fit from iRobitaille et al.l (|2007l ) for IRc9 is shown as the black line in the 
left panel. The dashed line shows the stellar photosphere corresponding to the central source of the 
best fitting model, as it would look in the absence of circumstellar dust (but including interstellar 
extinction). Right Panel — SE P for LVl, including our mid-IR FORCAST measurem ents, Spitzer 



IRAC fluxes (jMegeathl 120121 ) , and additional mid-IR ground-based measurements (jSmith et al 



2OO5I ). Photometry for LVl at 19.7 /xm was done using a 3.8" aperture. Upper limits (triangles) 



are 3-0" and based on RMS noise in the background at the position of the source. 
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perature throughout the central part of the nebula is ~70K and the emission is optically thin. The 
Ney- Allen Nebula is the warmest region in the nebula (besides BN/KL), with a color temperature 
in excess of Ri90K (see discu s sion i n Section [3. ip . These temperatures are significantly colder than 
those reported by iHavwardl (|1994| ) for both the warm (T ~ 300 K) and cold (T ~ 110 K) dust 
components. T his is to be expected, however, since our study utilizes longer wavelengths than 



HaywardI (jl994l ) and hence is more sensitive to colder dust. 



Discussion 



3.1. Ney-AUen Nebula 



The dominant feature of the Trapezium region is the Ney-Allen Nebula, which is bright in all 
the FORCAST filters except for 37 //m, where it is very difficult to distinguish from variations in 
the background emission (Fig. [3|) . The morphology of Ney-Allen is clearly wavelength-dependent 
in our FORCAST images: At 7.7 fim, it is an elongated source (approximately 4" x 7"), offset by 
2.2 arcsec from 9^ D Ori (see Fig[T]); in the images at 19 and 31 ^m, it is a crescent surrounding the 
position of 9^ D Ori with the apex pointing roughly toward a PA of 220°. The crescen t shape is also 
apparent at N and Q bands (jRobberto et al.ll2005l ). and 11.7 fim ([Smith et al.ll2005l ). In addition, 
the diameter of the crescent increases with wavelength, from approximately 5.5" at 19 /im to 8.2" 
at 31 fim, suggesting a radial temperature gradient indicative of heating by the central source, 
9^ D Ori. The same increase in size with wavelength (exceeding the increase in PSF size as a 
function of wavelength) is also evident in IRAC images of the Trapezium available from the Spitzer 
Legacy Archive. 



Robberto et al.l (120051 ) suggested that the Ney-Allen emission cou M be due to dust th at was 



evaporated from a massive circumstellar disk associated with 6^ D Ori. ISmith et al.l ( 20051 ) argue 
against this interpretation citing lack of evidence for such a disk (though it is possible that the 
disk is now completely dissipated) . Since the axis of symmetry of th e Ney-Allen cresce nt aligns 



very closely with the direction of proper motion for 9'^ D Ori (217°). ISmith et al.l (|2005l ^ instead 
postulated that the Ney-Allen emis sion is due to dust s wept up in the bow shock between the stellar 
wind of 0^ D Ori - a B0.5V star ( Mason et al.lll998l ) - and the ambient nebular material. Our 
SOFIA-FORCAST observations are consistent with this interpretation. Because 9^ D Ori is the 
only member of the Trapezium to show such a compressed shell, then, if this scenario is correct, it 
is probably in a relatively high density portion of the overall nebula. This suggests that 6^ D Ori 
is much closer to the background dense cloud tha n are the other rn embers of the Trapezium, and 
may even be plunging into the background PDR (ISmith et al.ll2005l ). 




Fig. 5. — Temperature map of the Orion Nebula from 19.7 and 37.1 /im data. White crosses mark 
the locations of the Trapezium stars {9^ Ori). 
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3.2. Proplyds 



Two of the compact sources detected in the 7.7 /um hnages (HST3 and LVl) are known 
proplyds — externally illuminated, low-mass young stellar objects whose circ umstellar disks are 
being photoabla ted by the strong EUV and FUV radiation from 9^ C Ori ( O'Dell et al 



1993 



Bally et alJll998l ). LVl is also visible in the 19.7 and 31.5 fim images, but is indistinguishable from 
the nebular background at 37.1 fiui. At 19.7 fim LVl appears slightly elongated in the NE-SW 
direction (Fig. [2]); the same orientation as the wind-wind bow shock (the stand-off shock between 
the proplyd evap orative flow and t he stellar wind fro m 6^ C Ori) that i s seen in both HST [O III] 
and Ha images (IBally et al.lll998l ) and at 11.7 //m (jSmith et al.ll2005l ). This suggests that dust 
from the proplyd circumstellar disk is entrained in the evaporative flow and then "piles up" at the 
wind-wind shock where it is heated by the UV radiation from 6^ C Ori and emits strongly in the 
mid-IR. Both HST3 and LVl are detected in 3.6 cm radio contin uum (jZapata et al.ll2004l ): only 
HST3 is detected in the COUP X-Ray survey (jCetman et ahlboosl ). 



We have combined our mid-IR flux meas urements with IR AC fluxes (jMegeathl l2012l ) , and 
additional ground-based mid-IR measurements (jSmith et al.ll2005l ) to generate the mid-IR SED for 
LVl shown in Fig. HI Interpretation of this SED is difficult, however, for two reasons: 



1. LVl is a known double (iGraham et al.ll2002l . and references therein), so mid-IR dust emission 
can arise from a circumstellar disk around either source and from the observed interaction 
zone between the two sources; 

2. mid-IR emission from the wind-wind shock cannot be spatially separated from the binary 
proplyd emission. 



Models of externally illuminated disks constructed by iRobberto et al.l (j2002l ) show that the 
mid-IR emission from 10 to 40 /im is dominated by emission from the disk interior (peaking at 
~ 40 nm) and emission from the superheated atmosphere (peaking at ~ 25 fj,m). These models 
tend to peak around 30 - 60 fim, which is not observed for LVl (the peak in the LVl SED appears 
to be around 20 /^m), suggesting that the mid-IR emission is dominated by a combination of the 
superheated disk atmosphere(s) of the two LVl components an d the dust piled up in the wind bow 



shock. LVl also has a modest 10 /im silicate emission fe ature (IShuping et al 



a strong component of the superheated disk atmosphere (IRobberto et al.ll2002l ) 



20061 ) ■ which is also 



3.3. Embedded Young Stellar Objects 



Roughly 30^^ north of the BN/KL region is the relatively isolated infrared source IRc 9 (|Wvnn- Williams &: Beck 
I974I : bownes et all Il98ll ). IRc 9 is thou ght to be an embed ded, luminous young stellar object 



(YSO) with a massive circumstellar disk ([Smith &: Ballvll2005l ): The near- to mid-IR SED shown 
in Fig. |4] is consistent with this interpretation. Model SEDs of young stellar objects were fit to 
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the SED of IRc 9 usi ng the Online Y S Q SE P Fitting tool (|Robitaille et alJ 120071 1 with best fit 
shown in Fig. [H The iRobitaille et alJ (|2007l ) models assume a single central stellar source with 
different combinations of axisymmetric circumstellar disks, infalling flattened envelopes and out- 
flow cavities. The fitting procedure interpolates the model fluxes to the apertures used in the 
measurem ents, scaling them to a given distance range of a source, which in this case is taken to be 
416ib6 pc (iKim et al.l 120081 ). and the line-of-sight e xtinction Ay. We adopt ed a range of Ay from 



2.0 to 10.0, based on the OMCl extinction map (jScandariato et al 



2011 



While these models 

yield estimates of physical parameters which are highly degenerate in general, they are fairly good 
at estimating the total luminosity of sources, which is our primary interest. 

The best fit SED models for IRc 9 have a to tal bolometric l umino sity of 110 to 215 Lq, 
which is somewhat higher than the value found by ISmith &: Ballyl (j2005l ) , but agrees with their 
conclusion that IRc 9 is the progenitor of an early-A (or late-B) ZAMS star with M = 3 — 4 Mq. 
The best fit models also indicate that the bulk of the emission longward of 10 //m is due to a large 
circumstellar envelope {R > 2700 AU), whereas the flux from 3 — 10 /xm is primarily from a much 
smaller circumstellar disk (R < 100 AU). T his suggests that the spatially resolved 8.8 to 18.3 fim 
emission observed by ISmith &: Ballvl (|2005l ) and attributed to a large circumstellar disk, may in 
fact be the inner regions {R < 1000 AU) of an envelope. 

In addition to IRc 9 we also detect three of the infrared sources in the OMCl-S region (IRSl, 
IRS2, IRSIO), also thought to be YSOs (Fig. [T]). At longer wavelengths, detection and flux mea- 
surement of the OMCl-S sources is not clear due to confusion with variations in the background 
emission. Limiting fluxes based on variations in the nebular background, however, suggest that 
20 - 40 /im emission from the OMCl-S sources may be s ignificant. We compared the 19, 31, and 
37 /um limiting fluxes to flux measurements from Spitzer (jMegeathll2012l ) and found that the SEDs 
from 4.5 through 31 fj,m are consistent with Class I or "flat spectrum" protostars, though it is 
not possible to rule out the possibility that they are highly-reddened Class II sources. The bright 
7.7 /im source SW of the BN/KL region coincides with HH202, a shock-ex cited region thought to 
be associated with an outflow or jet from OMCl-S IRSl (jSmith et al.ll2004l ). In contrast, HH202 is 
very weak at 6.6 ^m and not detected in the 19, 31, and 37 ^m images at all, strongly suggesting 
that the 7.7 ^m excess is due to PAH emission. 



3.4. Diffuse Emission 



The diffuse emission observed in our FORCAST images is generally a combination of emission 
from dust in both the background MIF and the H II region itself. The observed 19.7/37.1 /im 
color temperature values across the nebula (Fig. [5]) are consistent with cool dust in the back- 
ground PDR. The structures apparent in our 19.7 and 37.1 fim imag es are consistent with higher- 
resolution ground-based Q-band (20 /um) observations obtained bv iRobberto et al.l (j2005l ). The 
filamentary st ructure SW of the Trapeziu m coincides roughly with bright [N II] and Ha emis- 
sion features (lO'Dell &: Yusef-Zadehl l2000l ) and also with 8.4 and 1.5 GHz radio continuum emis- 
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sion (jPicker et al.ll2009l ). iRobberto et alj (|2005l ) concluded that the filaments in the mid-IR cor- 
respond to undulations or folds in the background PDR surface that cause limb-brightening and 
shadowing effects — similar to the Orion Bar but on a smaller scale. 

The diffuse emission drops significantly in the OMCIS region in all of our FORCAST images. 
The 19.7/37.1 color temperature drops significantly as well, indicating reduced UV heating in this 
region. We interpret this as a shadowing effect wherein the bright folds and undulations NE of 
OMCIS correspond to dense material that blocks the UV radiation field from 9^ C Ori. This is 
consistent with t he interpretation of OMCIS as a cold, dense region of gas and dust embedded within 
the H II region (jO'Dell et al.ll2009l ). We suggest that the formation of the embedded protostars 
in OMCIS may have been triggered by advancing shock waves when the Orion Nebula was first 
forming. 

The 37.1 fim image presented here is the highest spatial resolution image of the Orion Nebula 
at this wavelength to date. There are a number of features in the 37.1 /im image that have no 
counterparts at shorter wavelengths. Most prominent, is a ~60" arc of emission to the west of 



OMC IS (see Fig. [3]). This feature is not evident in na rrow-band emissio n (jO'Dell &: Yusef-Zadeh 
2000), but is present in the 1.5 GHz radio continuum (jPicker et al.ll2009l ). The arc coincides with 
a region of low color temperature as well (Fig. [5]). Interpretation of this arc is unclear, but it seems 
plausible that it is heated internally by the embedded OMCIS protostars. 

Most of the other structu res observed at 3 7 fim coincide with bright region s of [Nil] and Ha 
emission, seen both in natural ([Ballv et al. I ll998l l and extinction-corrected images (|0'Dell &: Yusef-Zadehl 
2OOOII. In particular, there is a very bright extended source just NW of the Trapezium (also noted 
by IWvnn- Williams et al.l (| 19841 ) in their low-resolution 30 ^m image) which corresponds to one of 
the brightest regions in the nebula in Ha and [N II] emission (jO'Dell Yusef-Zadehl I2OOOI ). Some 
of the diffuse structures at 37.1 ^m are also coincident with fluctuations in the color temperature 
(Fig. [5]). These two phenomena suggest that these are all perturbations in the surface of the back- 
ground PDR. One possible exception, however, is the emission WSW of 6^ C O ri, which coincide s 
with roughly 5 known proplyds, including the prominent sources LV4, 5, and 6 (iBally et al.lll998l ). 
We suggest that this is emission from dust that is entrained in the photoablation flows off the 
circumstellar disks around these sources and then heated by the UV radiation field of 6*^ C Ori. 
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